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ABSTRACT 

We propose a novel method for determining the ages of low-mass, pre-main sequence stellar systems using 
the apsidal motion of low-mass detached eclipsing binaries. The apsidal motion of a binary system with an 
eccentric orbit provides information regarding the interior structure constants of the individual stars. These 
constants are related to the normalized stellar interior density distribution and can be extracted from the predic- 
tions of stellar evolution models. We demonstrate that low-mass, pre-main sequence stars undergoing radiative 
core contraction display rapidly changing interior structure constants (greater than 5% per 10 Myr) that, when 
combined with observational determinations of the interior structure constants (with 5 - 10% precision), allow 
for a robust age estimate. This age estimate, unlike those based on surface quantities, is largely insensitive to 
the surface layer where effects of magnetic activity are likely to be most pronounced. On the main sequence, 
where age sensitivity is minimal, the interior structure constants provide a valuable test of the physics used 
in stellar structure models of low-mass stars. There are currently no known systems where this technique is 
applicable. Nevertheless, the emphasis on time domain astronomy with current missions, such as Kepler, 
and future missions, such as LSST, has the potential to discover systems where the proposed method will be 
observationally feasible. 

Subject headings: binaries: eclipsing — stars: evolution — stars: low-mass 
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1. INTRODUCTION 

In the study of the structure and evolution of low-mass stars, 
there are a variety of different methods capable of yielding 
reasonably accurate age estimates; for a thorough discussion 
of the different methods, and their strengths and weaknesses, 
see the review by Soderblom (2010). One such method uses 
detached eclipsing binaries (DEBs) to assign an age. DEBs 
are fantastic systems for studying stellar evolution. Obser- 
vations can provide precise masses and radii for the compo- 
nent stars that are nearly model independent (see reviews by 
Andersen 1991; Torres et al. 2010). Tight constraints on the 
stellar masses and radii allow for stringent tests of stellar evo- 
lution models. Furthermore, the age of a DEB system can be 
derived if stellar models can predict the radius of each star in 
the binary at a common age and with a single chemical com- 
position. 

Deriving an age estimate from a DEB is straightforward 
once precise masses and radii are extracted from the data. 
However, the reliability of the age estimate is contingent 
upon the accuracy of the stellar models. Recently, pre- 
main sequence (pre-MS) and MS models of low-mass stars 
(< O.8M ) have received substantial criticism for not accu- 
rately predicting the radii of stars in DEBs. As the number 
of DEBs with precisely measured masses and radii has in- 
creased, it has become clear that stellar models under-predict 
the radii of stars in DEBs by upward of 10% (see, for ex- 
ample, Mathieu et al. 2007; Jackson et al. 2009; Torres et al. 
2010; Feiden & Chaboyer 2012a). The discrepancies be- 
tween model and observed radii have been largely attributed 
to the effects of magnetic fields and magnetic activity (Ribas 



2006; Chabrier et al. 2007; Jackson et al. 2009; Morales et al. 
2010). When present, the discrepancy between observations 
and model predictions severely limits the use of stellar evolu- 
tion models to derive the age of individual DEB systems. 

We propose a novel method to use low-mass, pre-MS DEBs 
to estimate the ages of young stellar systems. Instead of com- 
paring individual stellar surface properties to stellar evolution 
models, we propose to use the dynamics of the DEB system. 
That is, comparing the observed rate of apsidal motion to stel- 
lar model predictions computed using the interior structure 
constants. This technique is less sensitive to the surface ef- 
fects of magnetic fields than are methods that only invoke the 
stellar radius or effective temperature. Our method has the 
potential to provide a more reliable age estimate. 

Below, we outline the model calculations (Section 2), in- 
cluding the computation of the interior structure constants 
(Section 3). Results pertaining to the time evolution of the 
interior structure constants as a function of stellar mass are 
given in Section 4. We conclude with a discussion of the use- 
fulness of this method in Section 5. 

2. MODELS 

2.1. Microphysics 

Model evolutionary tracks used in this study were 
computed with the Dartmouth Stellar Evolution Program 
(DSEP). 1 The physics incorporated in the models have been 
described extensively in the literature (Chaboyer & Kim 
1995; Chaboyer et al. 2001; Bjork & Chaboyer 2006; 

1 Available at http://stellar.dartmouth.edu/models/ 
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Dotteretal. 2007, 2008; Feiden et al. 2011), but we will 
provide a brief overview of the physics pertinent to the 
present study. 

Arguably the most critical component of low-mass stellar 
evolution models is the equation of state (EOS). For masses 
considered in this study DSEP uses theFreeEOS 2 inthe EOS4 
configuration. We selected the FreeEOS for three primary 
reasons. First, it includes non-ideal contributions to the EOS, 
such as Coulomb interactions and pressure ionization, that be- 
come important in the dense plasma of low-mass stars. Sec- 
ond, the FreeEOS calculates the EOS for hydrogen, helium, 
and eighteen heavier elements as opposed to only calculat- 
ing the EOS for hydrogen and helium. Finally, the FreeEOS 
may be called directly from within the stellar evolution code. 
This feature avoids the need to interpolate within EOS tables, 
thereby minimizing numerical errors. 

Conditions near the outer, optically thin layers of low-mass 
stars preclude the use of gray atmosphere approximations 
(Chabrier & Baraffe 2000, and references therein). There- 
fore, we use the phoenix ames-COND model atmospheres 
(Hauschildt et al. 1999a,b) to define the surface boundary 
conditions for our interior models. The atmosphere models 
are attached at the photosphere, 3 taken to be the point where 

T = r eff . 

The radiative opacities we adopt are the OPAL opacities 
above 10 4 K (Iglesias & Rogers 1996) in combination with 
the Ferguson et al. (2005) opacities below 10 4 K. For models 
that are not fully-convective, helium and heavy element diffu- 
sion are treated according to the formulation of Thoul et al. 
(1994). Fully-convective models are assumed to be com- 
pletely and homogeneously mixed because the convective 
timescale is considerably faster than the diffusion timescale 
(Michaud et al. 1984). 

2.2. Solar Calibration 

The primary input variables of stellar evolution models 
are defined relative to the Sun. These input variables in- 
clude the stellar mass, the initial mass fractions of helium (Y,) 
and heavy elements (Z,), and the convective mixing-length 
(ocmlt)- Therefore, we must first define what constitutes the 
Sun for our model setup. We require that a 1 Mq model accu- 
rately predict the solar radius, the solar luminosity, the radius 
to the base of the solar convection zone, and the solar photo- 
spheric (Z/X) at the solar age (4.57 Gyr; Bahcall et al. 2005) 
By iterating over different combinations of Yj, Zj, and (Xmlt 
we are able to converge upon a solution for the Sun. The final 
set of variables that satisfies the above criteria for the solar 
heavy element composition of Grevesse & Sauval (1998) was 
Yi = 0.27491, Z, = 0.01884, and oc M lt = 1-938. 

3. INTERIOR STRUCTURE CONSTANTS 

The distribution of mass within a star in a close binary sys- 
tem is influenced by the star's rotation and by tidal interaction 
with its companion. Imagine two stars, A and B. The rota- 
tion of star A and the tidal interaction of star B with star A 
distorts the shape of star A. Instead of remaining spherically 
symmetric, the equilibrium configuration of star A will be- 
come ellipsoidal. Subsequently, the gravitational potential of 

2 By Alan Irwin: http://freeeos.sourceforge.net 

3 The choice of where to attach the model atmosphere to the interior model 
is an important one. Our experience indicates that attaching the atmosphere 
at T = r c ff is reasonable for M > 0.2M . ; ; for lower masses it is necessary to 
attach the model atmosphere deeper into the interior. 



star A will also become ellipsoidal. The same can be said 
from the perspective of star B. 

If the binary orbit is elliptical, the distorted gravitational po- 
tential will cause the orbit to precess. This precession may be 
likened to the precession of Mercury's orbit about the Sun — 
although, Mercury's precession is due to general relativity. 
The precession of the binary orbit is known as apsidal mo- 
tion. 

The rate of apsidal motion (cb; measured in degrees per cy- 
cle), or the rate at which the orbit precesses, is governed by the 
shape of the gravitational potential, which may be deformed 
as discussed above. Therefore, cb depends on the properties 
of the stars and of the orbit. Explicitly, 
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In the above equation, Q.k is the mean orbital angular veloc- 
ity, Q.j, nij, and Rj are the rotational velocity, the mass, and 
the radius of the z'-th component in the binary, respectively. 
Additionally, A is the semi-major axis of the orbit, 
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with e being the eccentricity of the orbit. Finally, the last term 
in Equation (1), k 2 , is the weighted interior structure constant 
observed for the two binary stars. In general, 



T- C2.1/Q.1 +C2,2&2,2 

k 2 = — — . 

C2,l +C 2 .2 



(5) 



Here, k 2 , i and £2,2 are the interior structure constants for each 
star. 

Equations (l)-(5) are derived from a y'-th order solid har- 
monic expansion of the gravitational potential (see, for exam- 
ple, Kopal 1978). The interior structure constant for a given 
star, k 2 , is the second-order term from a more general set of 
expansion coefficients, kj. These second-order coefficients 
quantify the central concentration — or radial distribution — of 
mass within a star (Kopal 1978). Lower values of k 2 corre- 
spond to a higher level of central mass concentration. Point 
sources, for example, have k 2 = 0. 

Observationally, we can not solve for the individual k 2 val- 
ues. However, in the above equations, every variable is a di- 
rect observable, except for k 2 . The latter must be inferred 
from observational determinations of cb and the c 2 j coeffi- 
cients. Since the individual k 2 { values depend on the stellar 
density distribution, they can provide deep insight into the va- 
lidity of stellar evolution models. To bypass the restriction 
that only k 2 can be inferred from observations, we use k 2 val- 
ues from stellar evolution models in combination with the ob- 
served c 2 i coefficients to derive a theoretical k 2 . 

DSEP is equipped to calculate the general interior structure 
constants, kj, for a star at every evolutionary time step. This 
is achieved by solving Radau's equation after each model it- 
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eration. Following the formalism outlined by Kopal (1978), 

a ^ + ^W to + 1} + ^ to -i)=JU+i) (6) 

with /' e {2,3,4,...} being the order of the solid harmonic, 
a is the radius of an equipotential surface (a — r when the 
surfaces are spherically symmetric), and 



a dZj 
Ej da 



(7) 



In the above equation, e 7 is the stellar deviation from spheric- 
ity. We also introduced p(a), the density of the stellar plasma 
at radius a, and (p), the volume averaged density at each ra- 
dius, 

2. [ a n ( n >\r.n, 
„3 
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If we assume that tides and rotation do not significantly alter 
the shape of a star, we are permitted to use a spherically sym- 
metric model to compute the interior structure constant. Our 
code employs a 4 th -order Runge-Kutta integration scheme to 
obtain a particular solution of Radau's equation at the stellar 
photosphere. Interior structure constants are then directly re- 
lated to the particular solutions at the surface, T| j(R), through 



,/ + l-Tb-(*) 

2 [j + r\j(R)\ ' 



(9) 



But, is it reasonable to assume that the stellar mass is not sig- 
nificantly redistributed due to rotation and tides? 

The effect of rotation on the central mass concentration of 
stars with a total mass greater than 0.8M Q was investigated 
in several previous studies (Stothers 1974; Claret & Gimenez 
1993; Claret 1999). For stars less massive than O.8M , the 
effect of rotation should be negligible because they have a 
higher mean density compared to solar-type stars. To test 
this assumption, we used Chandrasekhar's analysis of slowly 
rotating polytropes to estimate the amount of oblateness — or 
deviation from spherical symmetry — rotating low-mass stars 
may be expected to have. 

Chandrasekhar (1933) derived an analytical expression for 
the stellar oblateness for slowly rotating polytropes. The 
oblateness was defined to be the relative difference between 
the equatorial radius and the polar radius, 



7 = 



r eq *"pole 



(10) 
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with r eq and r po i e being the equatorial and polar radius, re- 
spectively. Polytropes were considered slowly rotating when 
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where Q. is the stellar angular velocity, G is the gravitational 
constant, and p ( is the mass density in the stellar core. Does 
this criterion apply to real low-mass stars? If we assume a 
rotation period of 1.0 day and a core density of 10 g cm~ 3 
(realistic for pre-MS low-mass stars), then % w 10~ 3 . 

The result of Chandrasekhar's analysis was a relation be- 
tween the stellar oblateness and % for different values of the 
poly tropic index, «, 



f 5.79x for/i = 1.5 
J = I 9.82x forn = 2.0 
141.8X forn = 3.0 



(12) 



Stars with masses below 0.65 M© are best represented by a 
polytrope with 1.5 <n < 2.0. Assuming « = 2.0, a low-mass 
star will have 7 ~ 0.01. This treatment indicates that the 
effect of rotation on the sphericity of low-mass stars is a 1 % 
effect; rotation will not be addressed in this study. 

Assessing the influence of tides is difficult. We expect 
spherically symmetric models to provide accurate estimates 
of the mass distribution if -c ^ roc he> where ^ roc he is the 
Roche lobe radius. 

We also caution that the configuration of the binary must 
be considered. If the rotational axes are not aligned, as found 
with DI Her (Albrecht et al. 2009), then the validity of the as- 
sumptions required to derive of ki no longer hold. See Section 
5 for a further discussion. 

General relativistic distortion of the gravitational potential 
also plays a role in determining the rate of apsidal motion 
(Gimenez 1985). This contribution to the apsidal motion rate 
can be added to the classical contribution (i.e., Equation (1)) 



<i>tot = a>iV + a>GR, 



(13) 



where cb tot , and <x>gr are the total apsidal motion rate, 
classical apsidal motion rate, and the rate predicted from gen- 
eral relativity. The general relativistic contribution does not 
depend on the mass distribution of the stars. Thus, general 
relativity does not affect the theoretical derivation of the in- 
terior structure constants, but must be accounted for prior to 
comparing the theoretical and observational determinations of 

h. 

4. RESULTS 
4.1. Single Stars 

Individual, solar-metallicity evolutionary tracks were com- 
puted for several masses ranging from the fully-convective 
regime (0.25 M®) up to masses where stars have thin con- 
vective envelopes (0.75 M®). Near the boundary where the 
transition from a radiative core to a fully-convective interior 
is expected (~ 0.35 M®), a finer grid of mass tracks was gen- 
erated to allow for further exploration. The evolution of £2 
with age for each mass track is presented in Figure 1 . The full 
collection of mass tracks (including analysis routines) used 
hereafter have been made available on online 4 . 

Stellar models that develop a radiative core have a rapidly- 
changing interior structure constant between the age of 10 and 
100 Myr. This can be observed in Figure 1(a). As the convec- 
tion zone recedes, the central regions of the collapsing pre-MS 
star create a more centrally-concentrated mass profile, lower- 
ing the derived value of k%. The result is that the value of 
Ic2 decreases by about 5% - 10% every 10 Myr for masses 
above 0.45 M®. Masses below approximately 0.45 M un- 
dergo variations up to about 5%. 

This period of rapid contraction continues until a small 
convective core develops, producing a star with three energy 
transfer zones: a convective core, a radiative shell, and a con- 
vective outer envelope. The star settles onto the MS once the 
small convective core subsides with the equilibration of 3 He 
burning. In Figure 1, this process is manifested by the upward 
turn of the mass tracks near 100 Myr, followed by a flattening 
of &2 as the star enters the MS. 

After first developing a small radiative core on the pre-MS, 
stars with masses between 0.29 M® and 0.35 M® eventually 
maintain a fully-convective interior (see Figure 1(b)). These 

4 https://github.com/gfeiden/k2age/ 
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Figure 1. The time evolution of the interior structure constant, ki, for stars of various masses. Stars that develop a radiative core during the pre-main sequence 
exhibit rapidly decreasing kj values between roughly 10 and 100 Myr. (a) Total range of masses considered in this study shown in increments of 0.05 Mq. (b) A 
detailed view of the transition to fully-convective interiors. Only the 0.37 M, ., model does not ultimately become fully-convective. 
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Figure 2. The influence of various model properties on the predicted evolution of fc? for two different stellar masses. From left to right, the properties investigated 
are: scaled-solar metallicity in (a) and (d), helium mass fraction (Y) in (b) and (e), and artificial radius changes in (c) and (f). Models in the top panels ((a)-(c)) 
show a 0.55 Mq star while the bottom series ((d)-(f)) show a fully-convective 0.25 U . star. Note, the legend for panels (c) and (f) is split between the top and 
bottom panels due to space restrictions. The separate tracks indicated by both legends are presented for each mass. 



small radiative cores manifest themselves as small dips in the 
evolution of ki, just as for stars that maintain a radiative core 
on the MS. While the rate of change in ki is rapid during the 
core contraction, the relative change in ki with stellar age is 
small (< 5%). It is evident from Figure 1(b) that near the 
fully-convective transition, stars that eventually end up with 
fully-convective interiors exhibit a degeneracy in the age-A:2 
plane. There is no differentiating between their pre-MS core 
contraction and the eventual reduction of the radiative shell 
using hi alone. Below 0.29 M©, a radiative core does not 
develop on the pre-MS. 

Additional models were generated to investigate the effect 
of specific stellar properties on the predicted values of ki. 
In Figures 2(a)-(f) we illustrate the results of changing the 
scaled-solar stellar metallicity, helium abundance (Y), and the 
effects of artificially inflating and deflating the stellar radius. 
For this exercise, two masses were selected to study the ef- 
fects on the two broad categories of low-mass stars: fully- 
convective stars and stars with radiative cores. 

It is apparent from Figure 2 that changes to the chemical 
composition have the largest effect on the value of the interior 
structure constant on the pre-MS. Variations in metallicity of 
0.2 dex translate into a 5% difference in the calculated ki for a 
0.55 M s star at a given pre-MS age. Similarly, large variations 
in the helium abundance have the ability to produce changes 
in ki at the 5% level. 

Without some prior knowledge of the stellar composition, 
the effects of such variations may be confused with a differ- 
ence in stellar age. The effects on ki are lessened as mass de- 
creases, until variations nearly vanish in the fully-convective 
regime (panels (d) and (e) of Figure 2). Note that a change 
in Y pushes ki in the same direction in both mass regimes, 
whereas a change in metallicity produces a change in one di- 
rection for the radiative core case but the opposite direction in 
the fully-convective case. While this behavior does not lift the 
age-composition degeneracy entirely, it could prove a useful 
diagnostic in DEB systems whose components straddle the 
fully-convective boundary. 

We attempted to mimic the possible effects of magnetic 
fields on the structure of our models by computing models 
at solar-metallicity with ocmlt = 0.5, 1, and 3. The modi- 
fied mixing-length represents magnetic suppression of con- 
vection in the deep interior. Additionally, models were run 
where we artificially changed the surface pressure by ±20%. 
The altered surface pressure represents the possible influence 
of star spots on the stellar photosphere. The magnetic Dart- 
mouth models (Feiden & Chaboyer 2012b) were not used as 
they have yet to be evaluated for stars on the pre-MS. 

Figures 2(c) and 2(f) show how h varies with changes to 
OCmlt and the surface pressure. During the pre-MS contrac- 
tion, a 20% change in the surface pressure results in a ~ 2% 
change to the stellar radius. The accompanying change in ki 
was found to be 0.2% and 1.0%, for a positive and negative 
change to the surface pressure, respectively. In the 0.55 M s 
models, increasing OCmlt to 3 yielded a stellar radius that was 
2% smaller than our solar-calibrated model with ki variations 
under 1% throughout the star's pre-MS contraction. Decreas- 
ing (Xmlt, however, produced larger variations in ki. At an 
age of 60 Myr the model radii appeared inflated by 5% and 
15% for OCmlt =1.0 and 0.5, respectively. The correspond- 
ing changes in ki were, respectively, 4% and 12%. Panel (f) 
of Figure 2 indicates that the 0.25 M Q , fully-convective mod- 
els experience the greatest differences at the youngest age and 
that these differences diminish until the different tracks con- 
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Figure 3. The evolution of the theoretical ki for a binary having a 0.55M 
primary and a 0.25M.., secondary. The "observed" stellar radii are fixed at 
0.62Rq and 0.41i? Q for the primary and secondary, respectively (see Section 
4.2). The orbit was chosen to have an eccentricity e = 0.2. 

verge on the MS. These changes to OCmlt and surface pressure 
are for illustrative purposes only. It is not at all clear that sim- 
ply reducing OCmlt is a suitable approximation for the pres- 
ence of an interior magnetic field nor that altering the surface 
pressure is a good approximation for star spots at the surface. 

4.2. Binary Systems 

The rapid evolution of ki for a single star can provide an 
accurate age estimate for that star, but what about for a binary 
system? We stated in Section 3 that observations are only 
able to provide the weighted mean value of ki for two stars in 
a binary. To derive an age estimate for a binary system, we 
must find the theoretical k%. 

Temporal evolution of the theoretical ki values can be ob- 
tained by combining two ki mass tracks using Equation (5). 
Computation of the ci,i coefficients requires precise knowl- 
edge of the stellar masses, radii, and the orbital eccentricity 
(see Equation (2) and note that A becomes irrelevant). The an- 
gular velocity term in Equation (2) can either be measured or 
approximated using the orbital eccentricity, assuming pseudo- 
synchronization (Kopal 1978), 



Hence our reason for focusing on DEBs: they can yield pre- 
cise estimates of the stellar and orbital properties. 

The evolution of hi is simplest for an equal mass binary. 
In this case, both stars contribute equally to hi, meaning the 
ki track is exactly equal to the two individual ki tracks. The 
discussion from Section 4. 1 on single star tracks would then 
apply to the binary system. 

A binary with unequal mass components is not so simple. 
Two mass tracks are required — one computed for the mass 
of each star — and must be combined as a single track using 
Equation (5). How does this effect our ability to extract an 
age estimate? 

We have provided an example of a weighted ki evolution- 
ary track in Figure 3. The masses selected are those of the two 
stars presented in Figure 2, M\ = 0.55.1/ . and Mi = 0.25.1/ . . 
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For this example, we arbitrarily adopted an orbit with an 
eccentricity e = 0.2. Weighting of the theoretical ki val- 
ues is insensitive to the semi-major axis. The stellar radii 
(R\ = O.62/?0 and R2 = 0.41/?©) were selected from a so- 
lar metallicity mass track at an age of 40 Myr. To simulate 
the type of data available to an observer, we elected to fix the 
radius, and thus, fix the C2j coefficients (used to weight the 
value of kj) at each age in Figure 3. 

The rapid evolution of hi that is taking place in the O.55M 
star largely dominates the relatively slow hi evolution of the 
O.25M star. While the weighted kj value for the binary does 
not evolve as rapidly as for a single O.55M0 star, the weighted 
value still changes by about 5% every 10 Myr. This is com- 
parable to models of single stars and should not significantly 
hinder any age analysis. 

5. DISCUSSION 

DEBs provide an excellent laboratory for testing stellar 
structure and evolution theory in different mass and evolution- 
ary regimes. The mass-radius plane is the strictest test of stel- 
lar models because these two quantities are best constrained 
by the observations (Torres et al. 2010). Results from studies 
performing such comparisons have led to the consensus that 
standard stellar evolution models are unable to accurately re- 
produce the observed stellar radii for masses below ~ Q&Mr., 
(e.g., Torres et al. 2010; Feiden & Chaboyer 2012a). The 
model inaccuracies are particularly evident among pre-MS bi- 
naries (Mathieu et al. 2007; Jackson et al. 2009). 

Radius discrepancies of approximately 10 - 15% are rou- 
tinely quoted between pre-MS DEBs and models. This makes 
it difficult to derive an age with less than 50% uncertainty 
from stellar positions in the mass-radius plane. Magnetic ef- 
fects, particularly surface spots, are thought to belie the ob- 
served radius deviations. Canonical stellar evolution models 
are non-magnetic and are therefore unable to properly account 
for magnetic modifications to convective energy transport and 
for the presence of magnetic spots on the stellar photosphere. 

We therefore advocate the inclusion of the interior structure 
constant, fe, to overcome these age determination inaccura- 
cies whenever possible. While this study does not lead to us 
to conclude that individual kj_ values (and thus ki) are entirely 
insensitive to magnetic activity in low-mass stars, it is evi- 
dent that £2 has the potential to be a better diagnostic of the 
pre-MS evolutionary state than the stellar radius when surface 
magnetic activity is present. 

For instance, setting (Xmlt = 0.5 increases the radius and 
&2 of a single star by 15% and 12%, respectively. Fixing the 
radius to determine an age leads to an age that is upward of 
180% greater than if we assume a solar calibrated OCmlt- On 
the other hand, fixing ki leads to only a 20% greater age. The 
age errors are then compounded when we consider both stars 
in the binary. The decreased sensitivity of the individual h% 
values makes the mean ki a superior choice compared to sur- 
face quantities like the radius and effective temperature. 

The age precision returned from observational determina- 
tions of k2 is dependent on the precision with which the obser- 
vational ko and the system's metallicity are known. For exam- 
ple, given an equal mass binary with 0.5 Mq stars, knowing ki 
with 5% uncertainty and the metallicity to ±0.2 dex yields a 
pre-MS age with an uncertainty of approximately 33%. Con- 
straining the metallicity uncertainty to ±0. 1 dex improves this 
age uncertainty to 20%. Furthermore, to obtain an age with 
5% precision would require ki to be measured with near 1% 



precision for a metallicity known to within 0. 1 dex. 

5.1. Observational Considerations 

The results presented and discussed above show that it is 
possible to precisely derive the age of a binary system from 
measurements of apsidal motion. However, obtaining precise 
observations of apsidal motion and measuring ki with 5% pre- 
cision is a painstaking task. The binary system must meet 
several criteria and the data must be of high quality. 

Foremost is that the binary must be a double-lined, eclips- 
ing system. While this has been alluded to, we have yet to 
illuminate precisely why this is so. The reason for requir- 
ing a DEB stems from the need for extremely accurate and 
precise stellar and orbital properties. Equations (1) and (2) 
reveal that derivation of k% requires exquisite knowledge of 
the stellar mass ratio, the stellar radii, and the orbital proper- 
ties (eccentricity and semi-major axis). Only data from DEBs 
can provide these quantities in a (nearly) model-independent 
fashion (Andersen 1991; Torres et al. 2010). 

A meticulous examination of the binary light curve and ra- 
dial velocity curve is of the utmost importance. The light 
curve must be assembled from multi-epoch, time-series dif- 
ferential photometry and must provide nearly complete phase 
coverage. This latter feature is essential. Not only must the 
primary and secondary eclipses be captured, but also the be- 
havior of the light curve out of eclipse. Such data is becoming 
increasingly available with the latest generation of photomet- 
ric surveys (e.g., Kepler, C0R0T 5 , SuperWASP 6 ) and those to 
come (e.g., LSST 7 , BRITE 8 Constellation Mission). 

Similar to the photometry, a large number of high disper- 
sion, high signal-to-noise spectra are needed to construct a de- 
tailed radial velocity curve. Attempts must be made to provide 
adequate phase coverage (see Andersen 1991) and uncover 
deviations due to the Rossiter-McLaughlin effect (Rossiter 
1924; McLaughlin 1924) for reasons we shall discuss mo- 
mentarily. Torres (2012) points out that the spectroscopy is 
still the limiting factor of quality in DEB analyses. 

Only with the quality of data described above and within 
Andersen (1991) and Torres et al. (2010) can a truly adequate 
analysis of a DEB be performed. However, acquisition of data 
of such quality is very rewarding and allows for a rigorous ex- 
amination and characterization of the stellar system. It would 
permit the measurement of the stellar masses and radii with 
extreme precision (below 2%). This is imperative consider- 
ing that Equation (2) depends upon the fractional radii (R/A) 
to the 5 th power. Additionally, actual measurement of 0) re- 
quires careful monitoring of eclipse times of minimum. This 
can only be performed if one has a densely populated light 
curve. 

The data would also permit measurement of the binary ec- 
centricity and semi-major axis with high precision. These 
properties affect both Equation (2) and the general relativistic 
contribution discussed in Section 3 (Gimenez 1985). Recall 
that the contribution from the latter must be removed from the 
total apsidal motion rate to derive the classical contribution 
given in Equation (1). Detailed radial velocity curves not only 
provide accurate mass and eccentricity estimates, but may 
also be used to investigate the inclination of the system. The 

5 Convection, Rotation, & planetary Transits 

6 Wide Angle Search for Planets 

7 Large Synoptic Survey Telescope 

8 Bright Target Explorer 
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theory presented in Section 3 relies on the assumption that the 
rotational axes of the two stars be parallel to one another and 
perpendicular to the orbital plane. It is possible to evaluate 
this restriction by using the Rossiter-McLaughlin effect in a 
manner similar to that presented for DI Her (Albrecht et al. 
2009). A detailed radial velocity curve may, additionally, be- 
tray the presence of a third body. The presence of a tertiary 
may affect the binary orbit, altering the derived apsidal mo- 
tion. 

One added benefit is that lengthy observations may help to 
identify — and thus correct for — the impact of star spots and 
magnetic activity on the eclipse profiles. Star spots have the 
ability to distort the light curve, which can diminish the accu- 
racy of the derived stellar properties (Windmiller et al. 2010). 
Removing the effects of spots is critical to obtain not only 
precise but also accurate stellar properties. The easiest means 
of obtaining detailed time-series photometry is through space- 
based satellites, such as CoRoT and Kepler. However, long- 
term ground-based observational efforts are beginning to pro- 
duce apsidal motion detections (e.g., Zasche 2012), demon- 
strating that it is feasible to cany out the necessary observa- 
tions using ground-based telescopes. 

Finally, by acquiring a large number of quality spectra, it 
may be possible to extract the projected rotational velocities 
(vsin;) and a modest estimate of the chemical composition. 
Spectroscopic determinations of cool star metallicities is no- 
toriously complicated, but most pre-MS DEBs will likely re- 
side near or in a cluster from which metallicity estimates can 
be extracted using the higher mass stellar population. In the 
event there is no known association from which to draw a 
metallicity, techniques based on low- and medium-resolution 
spectra are encouraging (see, e.g., Rojas-Ayala et al. 2012). 
While the validity of such techniques along the pre-MS is un- 
clear, they provide a viable starting point and typically pro- 
duce metallicities with uncertainties below 0.2 dex, the limit 
we recommend. 

5.2. Limitations 

The usefulness of k% as an age estimator is limited to the 
pre-MS, in particular, during the evolutionary phase where 
the radiative core is rapidly contracting. This typically cor- 
responds to an age between 10 Myr and 100 Myr (see Fig- 
ure 1). At the 5 - 10 % measurement level, this technique is 
also restricted to binaries where one of the stars has a mass 
above ^0.40 M . This restriction ensures the radiative core 
contraction and relative change in ki for the more massive star 
is rapid enough to dominate the theoretical ki evolution. Re- 
ducing the observational uncertainty in ki below 5% enables 
a more accurate age estimate and would allow for DEBs with 
lower mass components to be reliably analyzed. 

Precise measures of apsidal motion and metallicity are chal- 
lenging to obtain, but are already feasible and should only 
improve over time. Instead of observational limitations, the 
greatest limiting factor for using ki as an age indicator is the 
circularization of the DEB orbit. Mutual tidal interactions will 
circularize binary orbits over time (Zahn 1977). Apsidal mo- 
tion requires an elliptical orbit. If an equal-mass binary is 
to maintain an elliptical orbit for the duration of its pre-MS 
contraction, the orbital period must be at least 2.4 days (Zahn 
1977). The probability of discovering an elliptical binary de- 
creases with time, but this provides an additional consistency 
check. The age suggested by £2 should not be significantly 
older than the orbital circularization timescale. 



Finally, the age estimate is only as accurate as the stellar 
models. Verifying that a given stellar model produces the 
proper mass distribution, hence ki, may at first seem rather 
unreasonable. However, at least one known system, with the 
possibility of a second (KIC 002856960; Lee et al. 2012), is 
capable of providing validation of the physics incorporated in 
low-mass stellar evolution models. Carteret al. (2011) have 
indicated that by the end of the nominal Kepler mission, they 
will know the interior structure constants of KOI- 126 B and C 
with about 1% precision. Interior structure constants known 
with this precision can place stringent constraints on the equa- 
tion of state of the stellar plasma (Feidenetal. 2011). The 
veracity of low-mass models, and therefore the validity of 
their predicted interior structure constants, may be assessed 
according to the results from KOI- 126 and similar systems. 
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